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ABSTRACT 

Chromospheric evaporation refers to dynamic mass motions in flare loops as a 
result of rapid energy deposition in the chromosphere. These have been observed 
as blueshifts in X-ray and extreme-ultraviolet (EUV) spectral lines corresponding 
to upward motions at a few tens to a few hundreds of km s“^. Past spectroscopic 
observations have also revealed a dominant stationary component, in addition 
to the blueshifted component, in emission lines formed at high temperatures 
(~10 MK). This is contradictory to evaporation models predicting predominant 
blueshifts in hot lines. The recently launched Interface Region Imaging Spec¬ 
trograph (IRIS) provides high resolution imaging and spectroscopic observations 
that focus on the chromosphere and transition region in the UV passband. Us¬ 
ing the new IRIS observations, combined with coordinated observations from the 
EUV Imaging Spectrometer, we study the chromospheric evaporation process 
from the upper chromosphere to corona during an XI.0 flare on 2014 March 29. 

We hnd evident evaporation signatures, characterized by Doppler shifts and line 
broadening, at two flare ribbons separating from each other, suggesting that chro¬ 
mospheric evaporation takes place in successively formed flaring loops throughout 
the flare. More importantly, we detect dominant blueshifts in the high tempera¬ 
ture Fe XXI line (~10 MK), in agreement with theoretical predictions. We also 
hnd that, in this hare, gentle evaporation occurs at some locations in the rise 
phase of the hare, while explosive evaporation is detected at some other loca¬ 
tions near the peak of the hare. There is a conversion from gentle to explosive 
evaporation as the hare evolves. 
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Introduction 


In the standard solar flare m odel, magnetic reconn ection takes place in the corona and 
releases a large amount of energy (jPriest fc Forbesll2002h . The released energy is transported 
downward to the lower atmosphere by non-thermal particles and/or thermal conduction. 
Most of the energy deposits in the dense chromosphere and heats up the local plasma. Due 
to a pressure imbalance, the heated chromosph eric plasma is driven up in t o the corona. This 


proce ss is called chromospheric evaporation flNeupertl Il968l: iHiravamal Il974l: lActon et ah 


19821 1. Based on momentum balance, there ar e also some chromospheric mater ials moving 
downward, called chromospheric condensation fjFisheiill987l: ICanfield et ahlll 99011 . The evap¬ 
orated hot plasma fills the newly formed flare loops that are then visible as post-flare loops 
in soft X-ray and extreme-ultraviolet (EUV) wavelengths. 

Chrornosphe ric evaporation can be detected in imaging and spectroscopic observations. 


Silva et al.l fll997(l reported that the footpoints of a flari ng loop in hard X -ray and soft X-ray 


images from Yohkoh moved toward the loop-top source. iLiu et al.l (120061) also found that the 


RHESSI hard X-ray sourc es in an M-c l ass fla re rose above the footpoints and merged into a 


loop-top source. Recently, iNitta et al.l (120121) observed an evaporation upflow motion of sev¬ 
eral hundreds of km s“^ in soft X-ray images from Hinode. There have been, however, only a 
few such imaging observations of chromospheric evaporation. In fact, our knowledge of chro¬ 
mospheric evaporation is mainly from spectroscopic observations. Mass flows related to evap¬ 
oration can be detected by Doppler shift measurements in chromospheric and coronal lines: 
upward motions generate blueshifts and downward motions give rise to redshifts in the lines. 


Antonucci et ah 

19821, 

1985; 

Antonucci & Dennis 

1983; 

Canfield et ah 

1987; 

Zarro & Lemen 

1988; 

Wiilser et ah 

1994; 

Dine et ah 

1996) using the Bent and Braee Crystal Spectrometer 


(BCS) on board the Solar Maximum Mission (SMM; lActon et ahl 119801) and Yohkoh/BCS 


-,-1 


have also been obtained in the 


hot XIX line (e.g.. 


Teriaca et ah 

2003, 

2006: 

Brosius & Phillins 

2004; 

Del Zanna et ah 


20061 ) using the Coronal Diagnostic Spectrometer (CDS; iHarrison et al.l 119951) on board the 
Solar and Heliospheric Observatory {SOHO). In the mean time, redshifts corresponding to 
downward velocities of tens of km s~^ have been measured in the chromospheric and tran¬ 
sition region (TR) lines fe.g.. IWulser et ah 1994 : Ding et ah 1995 : Czavkowska et ah 1999 


Teriaca et ah 2003 . 2006 : Brosius 2003 : Kamio et ah 20051: Del Zaun a et ah 2006 ). Recently, 


both blueshifts and redshifts have been detected in different emission lines using the EUV 
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Imag i ng Spectrometer (EIS: ICulhane et al. 


20091: IChen &: Dindl2010l: iLi &: Ding 


2011 


2007li on board Hin ode fe.g.. iMilligan &: Dennis 


Doschek et al.ll2013l) . 


According to the observed Doppler shifts, chromosp heric evaporation ca n b e divided into 
two ty pes: gentle evaporation and explosive evaporation ( Fisher et alll985a b c: Milligan et al 
2006al]bl) . When the hot coronal lines are blueshifted and the lines formed in the upper chro¬ 
mosphere and TR are redshifted, it is considered as explosive evaporation. When all the 
coronal, transition regional, and chromospheric lines show blueshifts, it is referred to as 
gentle evaporation. These two types of evaporation are also distinguished by the speed of 
evaporating plasma: gentle evaporation is characterized by tens of km s“^ (subsonic) up- 
flows, while explosive evaporation characterized by hundreds of km s“^ (supers onic) upflows 


( Zarro fc Lemenlll988l and references therein). Using hydrodynamic simulations, iFisher et al. 
( 1985a 3) reported that there may exist an energy flux threshold of about 10^° 


erg cm ^ s ^ 


(with a 20 keV energy cutoff ) between gentle and explosive evaporations. However, recently, 
Rubio da Costa et al.l (1201511 found that stochastically accelerated electrons wit h an energy 

could also cause explosive evaporation. In addition, iReep et al. 


flux <10^° erg cm ^ 


,-i 


(120151) claimed that the threshold between explosive and gentle evaporation is not a con¬ 
stant in the beam flux, but depends on the electron energy and duration of heating. Note 
that, besides being driven by no n-thermal electrons, both types of evaporation could also 


be driven by the r mal c o nduction (IZarro &: Lemenlll988l: lMilhganll2008l: iBattaglia et al.ll2009 


Gan et ahl Il99ll . Il992l : Ihiu et al.l 120091 ). In addition, both t ypes of 


transition between them are observed in certain flare e vents (iBrosius 


evapo ration could appear not only in large flares (e.g. , IMilligan et a. 


2010 ) but also in microflares ( Brosius fc Holman 2010 : Chen fc Dine! 201ol) . 


evapo ration and the 


20091). and explosive 


2006al: IVeronig et al. 


Although a large number of observational results are consistent with chromospheric 
evaporation models, there still exist some discrepancies between observations and models. 
Evaporation models predict that hot emission lines, such as Fe xxill, Fe xix, and Ca xix, 
should exhibit dominant blueshifts. In observations, however, these lines show a dominant 
stationary compo nent with only a weak blueshifted component indicating upflows of hun¬ 
dreds of km s“^ ( Antonucci et al. 1982 : Ding et al. 1996 : Milligan fc Dennis 2009 ). A few 
scenarios have been proposed to explain the observations: (1) the stationary plasma is from 


the top of flare loops where evaporated material is stopped ( Doschek &: Warren 20051): (2 ) 


the plasma is moving in the direction perpendicular to the line of sight (iFalewicz et al.ll2009l) : 
(3) the stationary component comes from the loop strands in which the evaporation process 
has ceased, while the weak blueshifted component originates from the newly formed strands 
with ongoing evaporation, and the observed spectrum i ncludes contr ibutions of both com¬ 
ponents in the same held of view along the line of sight ( Brosius 20131) : (4) the resolution of 
current instruments is not high enough to observe the evaporating plasma at an early time 
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flDoschek &: Warred l2005l) . Another problem concerns where the reversal from upflows to 
downflows takes place in explosive evapor ation. It is usually considered that this reversal 


occurs at the site of energy deposition fe.g.. Milligan &: Dennisll2nn9l:lLi Sz Dindl2011 ) . Evap - 


oration models predict that the reversal te mperature is below 1 MK flFisher et ah 


1985all . 


which was supported by some observations (IKamio et ahll2005l: iDel Zanna et al.ll200611 . How¬ 
ever, hi g;her reversal temperatures, such as 2 MK, or even 5 MK, were also found in obser¬ 
vations (iMilliganl 120081: iLi fc Ding||201111 . These observations pose challenges to evaporation 
models. 


The recently launched Interface Region Imaging Spectrograph (IRIS; iDe Pontieu et ah 


20141) provides high resolution imaging and spectroscopic data in the UV wavelength range. 
Due to the observing modes, it needs some fortune for the spectrograph to capture a solar 
flare. Fortunately, IRIS has captured more than 200 flares even including a few X-class 
flares as of June 2015. Furthermore, for some flares, coordinated observations with another 
spectrometer, EIS, which provides high resolution spectroscopic data in the EUV wavelength 
range, are available. Combining the UV and EUV spectroscopic data obtained by these two 
instruments, we study the chromospheric evaporation process from the chromosphere and 
TR to the corona during an X-class flare. In Section [21 we describe the flare observations; 
then we show our results of spectral analysis in Section [3l hnally we give a summary and 
discussions in Section 01 


2. Observations and Data Reduction 

The flare analyzed here is an XI.0 flare that occurred on 2014 March 29 in NOAA AR 
12017. This is the first X-class flare that IRIS has observed. Coordinated observations have 
also been obtained with EIS for this flare. Figure [H shows a summary of the observations. 
The top panel presents the GOES 1-8 A soft X-ray light curve, which peaks at 17:48 UT 
as indicated by the vertical solid line. The vertical dashed lines mark the times of multiple 
scans of IRIS and EIS over the flaring region. The bottom panels show the IRIS slit-jaw 
image in 1400 A, UV 1600 A image by the Atmospheric Imaging Assembly (AIA), and EIS 
He II 256 A intensity map. These images are taken at ~17:54 UT, several minutes after the 
flare peak time. We can clearly see the two flare ribbons in the images. Note that the field 
of view (FOV) of IRIS overlaps with that of EIS. 

IRIS obtains slit-jaw images and spectra in the near-UV (NUV, 2783-2834 A) and far- 
UV (FUV, 1332-1358 A and 1389-1407 A) passbands using a slit with a width of 0.33". 
The pixel scale is 0.167". IRIS observed an area of 167"x 174" in the slit-jaw 1400 A images 
with a cadence of 26 s for the XI.0 flare. The slit scanned over an area of 14"x 174" with 8 









































5 


steps (coarse steps, implying a spatial resolution of 2" across the slit) in each run that took 
75 s. IRIS observed the active region for several hours till 17:54 UT. Here we use the data 
in the last 25 minutes (from 17:30 to 17:54 UT), which cover the rise phase through part of 
the decay phase of the flare. During this time period, the IRIS slit scanned over the flaring 
region 19 times, which are indicated by the grey dashed lines in Figure [1] The IRIS level 
2 data are used in our study. The data have been processed with dark current subtraction, 
flat held correction, and geometric and wavelength corrections. For the spectra, the position 
of spectral lines suffers an orbital variation due to the temperature change of the detector 
and the change of spacecraft-Sun distance over the orbit. This line center variation was 
eliminated through the routine iris-orbitvar_corrJ2.pro in the standard IRIS software data 
package. The IRIS lines studied here are listed in Table [1] 

EIS acquires EUV spectra in the short-wavelength range of 170-210 A and long-wavelength 
range of 250-290 A with a spectral resolution of 0.0223 A pixeU^. The 1" slit of EIS was 
used to scan over the haring region with a step of 4". There are 11 steps in each run that 
took 2 minutes and 14 s, with a scanning area of 44" x 120" that overlaps with the IRIS FOV. 
EIS scanning repeated 13 times from 17:30 UT until the end of the coordinated observations 
with IRIS. The EIS data are reduced using the standard EIS software data package, with 
corrections for detector bias, dark current, hot pixels, and cosmic ray hits. The tilt of the 
slit is also corrected. In addition, we make a correction for the line position shift due to 
temperature variation of the spectrometer. The EIS lines used here are also listed in Table 

m 


To measure Doppler velocity, we use average line center to determine the reference 
wavelength of spectral lines. For the high temperature lines, e.g., Fe xxi and Fe xxill, which 
cannot be seen before the flare onset, we use the average line center over a loop-top region 
in the decay phase of the flare as the reference wavelength by assuming that the plasma in 
this region has reached a static state without signihcant mass flows well after the flare peak. 
With this method, we hnd the reference wavelength for the Fe xxi li ne to be 1354.126 A, 
which is consistent with the value 1354.106±0.023 A in lYoung et al.l ( 20151) . Alternatively, 


we can also use the relatively cool lines in the same spectral window, namely O l, to calculate 
the reference wavelength of Fe xxi. This assumes that the cool O l line is at rest before the 
onset of the flare; then we hnd the offset between the average line center of O l, observed in 
quiescent regions by IRIS, and the theoretical line center given in the CHIANTI database, 
and apply this same offset to the Fe xxi line. For comparison, this method yields the rest 
wavelength of Fe xxi to be 1354.124 A, which is very close to the value obtained in the 
former method. For the other warm and cool lines, i.e., Fe xvi, Fe xil. Si IV, He ll, C ll, 
and Mg ll, we use the average line center before the hare onset as the reference wavelength. 
The uncertainty in the measured velocities for all the lines is estimated to be less than 10 
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km s ^ considering the spectral resolntion and calibration errors of IRIS and EIS data. 

Co-alignment between IRIS and EIS images is made as follows. We first co-align the 
IRIS slit-jaw 1400 A image with the AIA 1600 A image, both of which show similar featnres 
(see Fignre[T]). Then we co-align the AIA 304 A image, in which the emission is mainly from 
the He ll line, with the EIS He ll intensity map. In addition, we also correct offsets (mainly 
in the Y-direction) between different spectral windows of EIS itself. Considering the spatial 
resolntion of IRIS and EIS, namely 2"x0.167" and 4"xl", respectively, the accnracy of onr 
co-alignment between IRIS and EIS is estimated to be abont 4"x4". 


3. Analysis and Results 
3.1. Analysis of the IRIS UV spectra 

IRIS provides high resolntion spectra in the FUV and NUV passbands, most of which 
were not well stndied for solar flares before. In the following, we first present the spectra 
along the slit (Fignre H]) and the line profiles at flare ribbon pixels (Fignres [SHS]) observed 
by IRIS to show its capability of diagnosing chromospheric evaporation. Then we present 
the evolntion of Doppler velocity along the slit (Fignre [7]) to show the temporal variation 
of chromospheric evaporation at the two flare ribbons. Here we analyze the hot Fe xxi line 
formed at ~10 MK, the Si IV line formed in the TR, and the Mg ll and C ll lines formed in 
the npper chromosphere. Note that the C ll line conld also be formed in the low TR. 

Fignre [2] shows the slit-jaw 1400 A images and the Fe xxi. Si iv, C ll, and Mg ll h&k 
spectra along the IRIS slit at different scan steps and different evolntion stages of the flare. 
The top panels give the slit-jaw image and the spectra at the seventh step with the position 
indicated by the white vertical line in the slit-jaw image jnst two minntes before the flare 
peak time. From the image, one can clearly see the two flare ribbons which are cnt across 
by the IRIS slit. Aronnd both ribbons, the hot Fe xxi line emission appears and shows 
some signihcant blneshifts, while the other cool lines show a strong red-wing enhancement. 
Note that in the brightest part of the ribbons, the Si IV spectra and the slit-jaw image are 
satnrated. A carefnl inspection also reveals existence of a large nnmber of narrow lines in 
all the spectral windows, especially near the flare ribbon. In the middle panels, we plot 
the slit-jaw image and the spectra at the sixth step of the slit at the flare peak time. One 
can see that there appears a slight separation between two flare ribbons compared with the 
earlier time. Near the sonth ribbon, the Fe xxi line shows very strong blneshifts of more 
than 200 km s“^, and the velocity decreases qnickly to nearly zero toward the north that 
is close to the flare loop-top region. However, the Fe xxi line intensity of the loop-top is 
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much greater than that near the ribbon. By comparison, at flare ribbons, the cool Si IV, 
C II, and Mg ll lines show a red-wing enhancement. The bottom panels show the slit-jaw 
image and the spectra at the third step in the decay phase of the flare. One can see that the 
two flare ribbons have a further separation in opposite directions. One can also And some 
short post-flare loops in the south around the slit as indicated by the white vertical line in 
the slit-jaw image. It is seen that the spectra in the loop-top region show some particular 
features: the hot Fe xxi line is almost unshifted though its intensity is greatly enhanced, 
while the cool lines exhibit some broadening, which may be due to superimposed line-of-sight 
velocities of cooling plasma. 

The general result revealed in Figure [2] is that the hot Fe xxi line shows significant 
blueshifts and the other cool lines show a red-wing enhancement at the flare ribbons. This 
can be considered as clear evidence of chromospheric evaporation, in particular, explosive 
evaporation. To display the evaporation signatures during the flare more clearly, we plot 
the line profiles at four ribbon pixels (named as NO, SI, S2, and S3) marked by the red bars 
in the slit-jaw images (Figure |2]), and measure the Doppler velocities. When measuring the 
Doppler velocity, we adopt different methods for different spectral lines. For the optically 
thin Fe xxi and Si iv lines, we use Gaussian function to fit the observed line profiles. Since 
there are several narrow and cool lines (shown in the top left panel of Figure [3], ‘u’ means 
unidentified) blending with the Fe xxi line in general, we adopt a multi-Gaussian function 
with nine components representing nine lines and assume that the lines from the same ion 
having the same Doppler velocity and width in order to fit the Fe xxi profiles. However, a 
double Gaussian function is applied to the line profiles that exhibit only G l blending with 
Fe XXI, such as the ones near the flare loop top. Most of the Si IV line profiles show a 
good single Gaussian shape, and we adopt a single Gaussian function to fit the Si iv line. 
However, there are some Si IV line profiles showing an asymmetry, mainly at the flare ribbon, 
so we also use a double Gaussian function to fit those line profiles and compare with the 
fitting results by a single Gaussian function. It shows that the velocity derived from a double 
Gaussian fitting (the redshifted component) is a little larger than the velocity from a single 
Gaussian fitting, but it does not change the sign of the velocity. The G ll and Mg ll lines are 
optically thick, which are formed in the chromosphere through a complex radiative transfer 
process. The core (ha or ka component) of the Mg ll h&k resonance lines is formed in the 
upper chromosphere just below the TR, whereas the wings (h 2 and k 2 ) are formed lower. In 
the quiet Sun (and pre-flare), the Mg ll resonance lines usually exhibit a central reversal at 
the core. In this flare, the central reversal is not present in most of the observed Mg ll h&k 
resonance lines as well as the G ll lines at the flaring locations. The line profiles exhibit a 
single emission peak, which appears to be redshifted with respect to the quiescent line center. 
Here we adopt a bisector or moment analysis method to estimate the line center shift of G ll 


and Mg ll. Because of a similarity between the two Mg ll (or C ll) resonance lines, we only 
present the results for one of the lines. The observed line profiles and the htting results at 
the selected four ribbon pixels are shown in Figures [3HH1 

In Figure [3l we plot the line profiles of Fe xxi, Si IV, C ll, and Mg ll k, together with 
the fitting curves or net emissions, at the north ribbon pixel NO that brightened in the rise 
phase of the flare. It is seen that the fitting curve of the Fe xxi line matches the observed 
one very well. The Fe xxi line exhibits a significant blueshift with a Doppler velocity of 117 
km s“^. The Si iv line, however, shows a red asymmetry. When using a single Gaussian 
function to £t the line prohle, we get a redshift with a velocity of 15 km s“^. If using 
a double Gaussian htting, we can obtain a larger redshift velocity of 55 km s“^ with an 
intensity ratio of 0.6. Similarly, the G ll line shows a redshift with a velocity of ~36 km s“^ 
by a bisector method (at ~40% of the peak intensity) or 37 km s“^ by the moment analysis 
on the hare net emission (subtracting the pre-hare emission). The Mg ll k line also displays 
a red asymmetry, seen from the moment analysis on its pre-hare (-1-5 km s“^), during-hare 
(-1-14 km s“^), and hare-net (-1-30 km s“^) emissions, respectively. 

Figure m shows the line prohles and htting curves or net emissions at the south ribbon 
pixel SI that brightened also in the rise phase of the hare. The intensity at this pixel is 
much stronger than that at pixel NO. Specifically, the intensity of Fe xxi at SI is three 
times the one at NO, and the Si IV and G ll lines seem saturated at SI (here we analyze 
the unsaturated Si iv and G ll line prohles at the closest ribbon pixel nearby SI instead). 
In addition, the Mg ll lines show a stronger emission at SI than at NO. At SI, we obtain 
a blueshift with a velocity of 64 km s“^ in the hot Fe xxi line, and some redshifts or a red 
asymmetry with velocities of tens of km s“^ in the cool Si iv, G ll, and Mg ll lines. 

Figure [5] presents the line prohles and htting curves or net emissions at the south ribbon 
pixel S2 that brightened at the hare peak time. At this pixel, the Si iv and G ll lines are also 
saturated, and similarly, we analyze the unsaturated prohles nearby instead. The intensity 
of Fe XXI at S2 is not as strong as the one at SI, but the Doppler velocity at S2 reaches 
as high as 214 km s“^, much larger than that at SI. In addition, the cool lines also exhibit 
redshifts with velocities of 20-45 km s“^. The Mg ll line also show a broad line prohle. All 
of these indicate that the evaporation is very strong at this pixel in the impulsive phase. 

In Figure [6l we display the line prohles and htting curves or net emissions at the south 
ribbon pixel S3 that brightened in the decay phase of the hare. It is seen that the intensities 
of all the lines at S3 have signihcantly decreased compared with the line intensities at the 
previous pixels in the rise and impulsive phases. The Fe xxi line also shows a decreased 
blueshift velocity of 30 km s“^. The other cool lines still show some red asymmetry. 
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The high temperature Fe xxi line, however, is symmetric and well htted by a single 
Gaussian function with a blueshift velocity. This suggests that, on the flare ribbon, the Fe 
XXI line is blueshifted as a whole. In particular, at S2 in the impulsive phase, the line shows 
a blueshift with a strong velocity of 214 km s“^. This is different from previous results, 
which only show a weak blueshi f ted component superimp osed on a dominant stationary 
component fe.g.. lDing et al.lll996l: iMilligan fc DennisI 1200911 . Our results are consistent with 
the prediction by theoretical models that hot plasma is massively driven upward in an im¬ 
pulsively heated flare loop. It is likely that the high spatial resolution of IRIS allows us to 
resolve the locations where energy deposition and strong dynamics take place in the atmo¬ 
sphere, whereas previous observations with low spatial resolutions may reflect a mixture of 
unresolved multiple loops with both blueshifted plasma and stationary plasma. 


Besides the line prohles at ribbon pixels, we also analyze the spectra at all the pixels 
along the slit. Figure [7] shows the time sequence of Doppler velocity, as well as non-thermal 
velocity derived from the line width, of the Fe xxi and Si iv lines, for pixels along the slit 
between the two green bars indicated in the bottom left panel of Figure |2l The map covers 
a spatial scale of 30" along the slit, which lies across the two ribbons separating from each 
other. Since the Fe xxi line emission is weak in many pixels especially at the north ribbon, 
to get a better signal-to-noise ratio, we rebin the data to sum up spectra in six IRIS pixels 
along the slit. We have conhrmed that rebinning the data does not affect signihcantly the 
measured Doppler shift distribution along the slit. Now the pixel size is 0.33" by 1" in Figure 
[71 Also note that, the black data points in the Fe xxi velocity panel represent the strong 
Doppler velocities that are far beyond the colorbar range (±150 km s“^). Some of these 
velocities even exceed 250 km s“^. In addition, in Figure [7] we only show the time sequence 
of velocity along the slit for the third step. For the other steps, either the observed data 
contain bad points (with all the intensity counts being a negative value of -200), or the Fe 
XXI emission is too weak to be identihed at the north ribbon. 


More interestingly, we hnd an evident pattern of ribbon separation in the velocity maps, 
as shown in Figure [71 The separation of two flare ribbons is an indication of progressive 
reconnection that forms flare loops successively and deposits energy at the feet of these 
newly formed loops. In this event, the Fe xxi velocity map shows a pair of blueshift fronts 
separating from each other at an apparent speed of about 20 km s“^, the same as the ribbon 
separation speed. These blueshifts, corresponding to upflow velocities in the range of 50- 
300 km s“^, are located at the bright fronts of flare ribbons, as indicated by the intensity 
contours of the Mg ll h (or Si iv) line. The Fe xxi velocity map also reveals that, on the 
ribbons, particularly the south ribbon, the region with signihcant blueshifts extends over a 
spatial scale of 2-3 arcseconds along the slit. At the same location, the blueshift diminishes 
in several minutes. In the south ribbon, after the bright front passes, the blueshift turns into 
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a redshift (with a velocity of 3-15 km s“^). Most of the redshift velocities are within the 
measurement uncertainty. It is likely that some of the redshift signatures are real, though 
we are not certain what might be the reason for the redshift in the hot Fe xxi line. 


Similar evolution patterns of the dynamics on the flare ribbon fronts are also observed 
in other velocity maps. In particular, signihcant redshifts, with a velocity of a few tens of 
km s“^, in the cool Si iv line are observed at locations of strong ribbon emission. (Note 
that the persistent weak blueshift below the north ribbon may not be related to the flare.) 
At these places, the maximum line broadening in both the Fe xxi and Si IV lines is also 
observed to be of the order of 100 km s“^ in terms of non-thermal velocity. These large non- 
thermal velocities on the flare ribbons may imply multipl e velocity comp onents in unresolved 
strands when energy deposition takes place. In fact, iMilliganI (120111) has found a good 
correlation between non-thermal velocity and Doppler velocity in the Fe xv and Fe xvi lines 
at loop footpoints, and the author interpreted the excess broadening as a superposition of 
flows (turbulence), presumably caused by chromospheric evaporation. The dynamic features 
diminish between the two ribbons. 


These results provide clear evidence of chromospheric evaporation at the sites where 
impulsive energy deposition takes place. The time scale of the evaporation signature (partic¬ 
ularly the blueshift) at one site is estimated to be 1-2 minutes (might be constrained by the 
observing cadence), and along the slit, the spatial scale of signihcant blueshift signatures is 
2-3 arcseconds. The presence of coincident blueshift in the hot line and redshift in the cool 
lines suggests the regime of explosive evaporation. This evaporation process proceeds well 
into the decay phase of th e hare, indicative of continuous energy rel ease/deposition from the 
start to the decay phase ( McTiernan et ah 1993 : Jiang et ah 20061) . Evaporation occurs in 
diherent loops formed successively by reconnection, therefore, it exhibits an apparent motion 
outward with a similar velocity to the separation velocity of hare ribbons. 


3.2. Analysis of the EIS EUV spectra 

Aside from the Fe xxi line (~10 MK), IRIS observations are mostly made in chromo¬ 
spheric and transition region lines with characteristic temperatures up to 0.1 MK. In com¬ 
parison, EIS provides EUV spectra that focus on diagnosing properties of transition region 
and corona plasma of temperatures from 0.05 MK to dozens of MK. To have a comprehensive 
picture of atmospheric dynamics from the chromosphere to the corona, we also select four 
emission lines from EIS, listed in Table [H to study the chromospheric evaporation process 
in the XI.0 hare. These lines. He ll, Fe xil, Fe xvi, and Fe xxill, form at 0.05, 1.3, 2.5, and 
13 MK, respectively. Figure [H] shows the intensity map in these four lines. No emission is 
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observed in the very hot Fe xxill line during the early phase of the flare. At 17:46 UT, the 
start of the major GOES soft X-ray peak, this line shows signihcant emission in the flaring 
region. The intensity maps of the other cooler lines also exhibit two flare ribbons, with the 
south ribbon much stronger than the north ribbon. In the intensity map of the coolest He 
II line, in particular, one can notice separation of the two ribbons, which is consistent with 
IRIS observations. 


In Figure ini we present the Doppler velocity map of the four EIS lines constructed from 
multiple scans. Here we only show the results for the pixels that have enough signal-to-noise 
ratio (mainly for the Fe xxill and Fe xvi lines). All the velocities are derived by htting the 
observed line prohles to a single Gaussian function. From the Doppler velocity map, we can 
see that the hot Fe xxill line shows strong blueshifts at the flare ribbons. The relatively 
cooler lines of Fe xvi, Fe xil, and He ll show a variety of blueshifts before the flare peak time. 
Note that some strong buleshifts in these cooler lines may be related to a hlament eruption 
(for details see iKleint et al.l 1201511 during the early time of the flare. From 17:46 UT, two 
minutes before the flare peak, on the flare ribbons, the cooler lines become mostly redshifted, 
while the hot Fe xxill line stays blueshifted. These manifest explosive evaporation, which 
conhrms the IRIS results. Note that, in the velocity map of the hot Fe xxill line, there also 
appear some patches with redshifts near the south ribbon, which might be caused by the 
uncertainty in velocity measurements. 


We note that almost all of the Fe xvi, Fe xil, and He ll line prohles show a single peak, 
and can be well htted with a single Gaussian function (see bottom panels in Figure [TII]) . 
However, some of the Fe xxill line prohles exhibit a weak blue asymmetry or even double 
peaks (mainly at the hare ribbon), though the majority still show a single peak of Gaussian 
shape (especially near the loop-top region). These diherent types of Fe xxill line prohles 
observed during the hare are presented in the top panels of Figure [TUI In the following, we 
alternatively adopt a double Gaussian function to ht the Fe xxill line prohles that display 
double peaks or weak blue asymmetry, as shown in the top left and middle panels of Figure 

m 


Figure ITT] shows the evolution of Doppler velocities along the slit in the four EIS lines. 
The EIS slit position (its eighth step) is nearly the same as the IRIS slit position (its third 
step) shown in Figure [71 Due to the diherent spatial resolutions as well as diherent observing 
times, there may be a relative shift of several arcseconds in both X and Y directions between 
the EIS slit and IRIS slit. In addition, the time scales of EIS and IRIS scans are diherent: 
EIS takes a longer time than IRIS for each scan, and EIS observations persist hve minutes 
longer before the end of the hare. These maps also show an apparent separation motion 
between two ribbons exhibiting Doppler shifts. However, the pattern is not as evident as the 





12 


IRIS Fe XXI measurement in Figure [3, which is likely due to the low spatial and temporal 
resolutions of EIS. Note that, the Fe xxill Doppler velocity map is derived from a double 
Gaussian fitting to the line prohles displaying double peaks (within the red rectangular box 
in the panel) or a week blue asymmetry, and the velocities represent the shifts of blueshifted 
components, which can exceed 200 km s“^. The mix of the blueshifted and stationary 
components in the EIS Ee xxill line, most likely due to the lower spatial resolution of EIS, 
explains why the blueshift signature is not as prominent as demonstrated in the IRIS Ee xxi 
line. 


3.3. Doppler Velocity versus Temperature in Multiple UV and EUV Lines 


Combining the multi-temperature UV and EUV lines from IRIS and EIS, we can study 
Doppler shift signatures in cool lines and hot lines to distinguish different types of chromo¬ 
spheric evaporation at different locations during the evolution of the XI.0 flare. We select 
four time-space locations marked in Figure [TTl For each location, the spectra is from an 
area of 4"x4". Three locations (1, 2, and 3) are on the flare ribbons and brightened before 
the peak of the flare. These places show prominent evaporation signatures. At the fourth 
location, evaporation has stopped after the peak of the flare. Doppler velocities in different 
lines are plotted as a function of temperature in Figure [T21 It is seen that, at location 1, 
all the spectral lines with different temperatures are blueshifted, which can be considered 
as gentle evaporation. At locations 2 and 3, the relatively hot lines are blueshifted, while 
the cool lines are redshifted, which is regarded as explosive evaporation. Note that, at loca¬ 
tions 2 and 3, the velocity in the Fe xxill line marked by a vertical bar is derived from the 
blueshifted component, when a double Gaussian htting is applied. The plots also suggest 
that, in the event of explosive evaporation, the reversal temperature from upflows to down¬ 
flows is likely below 1 MK at location 2 and close to 1 MK at location 3, although there 
is a large uncertainty due to lack of spectral lines coveri ng the temperatur e range from 0. 1 
to 1 MK. These results are similar to previous studies flKamio et al .1120051: lMilliganll2008h . 
Different flow reversal temperatures would imply that the energy deposition height, rate, or 
duration may be different at these two locations (iLiu et al.ll2009l) . 


These four locations are at different positions along the slit (location 4 overlaps partly 
with location 1). Their Doppler velocities plotted in Figure [12] are for different times at 
different stages of the flare. For locations 1, 2, and 3, chromospheric evaporation is in 
progress. Evaporation at location 1 is gentle evaporation which begins at the start of the 
major peak of the flare. Evaporation at locations 2 and 3 is explosive evaporation which 
occurs about two minutes later, close to the peak time of the flare. Therefore, in the flare. 
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different locations can show different types of evaporation, and the evaporation can vary 
from a gentle one to an explosive one as the flare evolves. These results indicate that energy 
deposition, reflected by the magnitude and types of chromospheric evaporation, varies with 
time and space during a flare. In more details, there are two possible reasons for the time- 
and space-dependent energy deposition: (1) the injected energy flux of the electron beam 
or conduction front that deposits the energy is time- and space-dependent; and (2) the 
magnetic and density structure of the atmosphere, which determines heating and dynamics 
of the lower atmosphere, may change from place to place. 


4. Summary and Discussions 

We have studied the chromospheric evaporation process during an XI.0 flare combining 
the IRIS UV and EIS EUV spectroscopic data. From the IRIS data with high temporal 
and spatial resolution and a good time coverage of the flare, we have detected dominant 
blueshifts in the high temperature Fe xxi line (~10 MK) at the flare ribbons. The Fe xxi 
line prohle at most locations on flare ribbons is symmetric and very well htted by a single 
Gaussian function, and the entire line is blueshifted. The dominant blueshifts throughout 
the flare have rarely been observed by other spectrographs. Therefore, our result provides 
an observational support to theoretical models predicting that high temperature lines should 
be entirely or mostly blueshifted. 

We have constructed Doppler velocity maps derived from multiple lines obtained by 
IRIS and EIS slits at multiple times during the rise, peak, and decay of the flare. We have 
found that blueshifts in the high temperature Fe xxi line and redshifts in other cooler lines 
from IRIS are located at the bright fronts of two flare ribbons, and spread out in the same 
manner as the separation of two ribbons. This evaporation process is seen to proceed into 
the decay phase of the flare. In addition, line broadening is also observed at flare ribbons. 
These results conhrm that dramatic dynamic processes occur in the flaring atmosphere when 
energy deposition takes place. The pattern of ribbon separation, together with the dynamic 
features, reflects continuous magnetic reconnection and energy release. Similar evolution 
patterns are observed in velocity maps constructed with the EIS EUV spectra, though with 
a much lower spatial (and temporal) resolution. 

By analyzing multiple lines observed by IRIS and EIS, we also examine Doppler shifts 
in cool lines and hot lines to study different types of chromospheric evaporation. In this 
flare, the evaporation process shows temporal and spatial variations: gentle evaporation is 
observed at one location during the rise of the flare, and explosive evaporation is observed 
at some other locations several minutes later, close to the peak of the flare. The conversion 
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from gentle to explosive evaporation dnring the XI.0 flare may be related to non-thermal 
electron beams. We have examined the X-ray light curves from RHESSI, which show evident 
non-thermal hard X-ray emission up to 100 keV. The non-thermal hard X-ray light curve 
peaks at 17:46 UT, close to the scanning time at locations 2 and 3, both of which show 
explosive evaporation. By comparison, location 1 shows gentle evaporation that occurs two 
minutes earlier when high energy hard X-ray emission has not been enhanced. 


Using high resolution spectra from IRIS, we find dominant blueshifts in the high tem¬ 
perature Fe XXI line. The evaporation model predicts that the high temperature lines should 
be entirely or mostly blueshifted. However, there have been very few spectral observations 
of such kind to support the model prediction. In the past, Doppler shift measurements were 
obtained in high temperature lines such as Ca xix by BCS and Fe xix by CDS, respectively, 
with low spatial resolution. These measurements showed a weak blueshifted component 
superimposed on a dominant stationary component. These earlier observations could not 
reveal a dominant blueshift component, perhaps for two reasons. First, dominant blueshifts 
in hig h temperature lines usually appear in the early evaporation phase (iDoschek fc Warren 


20051 1 when the high temperature emissions are very weak. Secondly, these dynamic features 


may occur in a very small region that could not be resolved with past instrument capabili¬ 
ties, and the observed spectra included contributions from nearby regions. Recently, the Fe 
XXIII and Fe xxiv lines with high spatial resolution from FIS were used to diagnose mass 


component ( 

Watanabe et ah 2010: Young et ah 2013) or that these lines are blueshifted as 

a whole ( 
example. 

BrosiusI 

2013 

). Nevertheless, such observations were too few to be conclusive. For 

Brosius 

(2013) found the above phenomenon appears at one location and lasts for 


only about two minutes. Here, we use the Fe xxi line from IRIS to study the evaporation 
pr ocess. We fad that the line is blueshifted as a whole (also reported in different fla res 
by Tian et ah 2014b . 2015 : Young et ah 2015; Polito et ah 2015 : Graham fc Cauzzi 2015 re¬ 
cently) in the early evaporation phase in multiple loops continuously formed during the flare, 
thanks to the high temporal and spatial resolution of IRIS, as well as its high sensitivity in 
detecting high temperature emissions in the early phase of flare heating. 


From the high resolution spectra by IRIS, we have observed a dynamically evolving 
evaporation pattern characterized by blueshifts in high temperature (~10 MK) lines, which 
are co-spatial and coincident with separating flare ribbons. Flare ribbon separation has been 
observed for many years in Ha, UV, EUV, as well as hard X-ray images. However, previous 
studies have rarely reported separation of flare ribbons accompanied by coincident significant 
blueshifts, in particular, in high temperature lines. Flare atmospheric dynamics have been 
extensively studied; most previous studies have reported time evolution of chromospheric 
evaporation at some locations or spatial distribution of Doppler velocities at some given 
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times (e.g.. ICzavkowska et al. 


Brosius fc Holman 


201 


I 


1999 


Brosius 


BrosiusI I 2 OO 3 I: iBrosius &: PhillipsI l2004 iBrosiusI 12009 


20131) ■ However, as limited by the time cadence and spatial 


resolution of the spectrometers, these past studies have not reported the spatio-temporal 
variation of the velocity pattern as shown in this work. Our result thus provides novel ev¬ 
idence of flare dynamics, which, together with what have been found in previous studies, 
can be explained in terms of the standard flare model; magnetic reconnection occurs suc¬ 
cessively at growing altitudes forming new and longer flare loops, and separating ribbons 
outline the feet of these successively formed loops, along which energy flux is deposited to 
heat the lower atmosphere and generate chromospheric evaporation at the feet; this process 
results in apparent separation of the velocity ribbons seen in the Dopplergrams. Note that 
reconnection proceeds continuously into the decay phase of the flare. 


We also note that the time and spatial scales of the strong blueshifts (in the hot Fe xxi 
line) and redshifts (in the cool Si iv line) at the ribbon are 1-2 minutes and 2-3 arcseconds, 
respectively. In general, these scales are a result of convolution of several factors: (1) the 
spread speed of flare ribbon (related to reconnection rate): the faster the spread motion 
is, the wider the blueshift band would be; (2) the duration of energy deposition at the loop 
footpoint, immediately after energy release by reconnection; (3) the hydrodynamic time scale 
within the flare loop; and (4) temporal and spatial resolutions of observations. Studying the 
ribbon dynamic features as revealed by emission features at different layers would help us to 
understand the flare heating source, rate, and duration, and further provide constraints to 
loop heating models and numerical simulations of flares. Note that the temporal and spatial 
scales of blueshifts/redshifts might be different in different flares. 


Besides high temperature lines, IRIS provides important chromospheric and TR lines 
that can help us diagnose the lower atmospheric condition during chromospheric evaporation. 
These lines include the Mg ll h&k, C ll, and Si IV resonance lines. The Mg ll and C ll lines 
are optically thick, and formed in the chromosphere through a complex radiative transfer 
process. Interpretation of these optically-thick lines requires sophisticated methods such as 
through radiative transfer modeling. In the present study, we adopt a bisector or moment 
analysis method on the flare net emission to estimate the plasma motion in the chromosphere. 
Our results provide an estimate of the velocity held in the relatively low layers (i.e., the upper 
chromosphere), in comparison with the velocities in the transition region and corona deduced 
with optically-thin lines. 


The Mg II and C ll line prohles become pure emission in the core from a central re - 
versal during the hare, which may suggest a greatly reduced opacity flTian et al.N2014al) . 
The opacity reduction might be caused by energy deposition in the line formation layer. 
Furthermore, a sharp change in the spectra in these lines (see Figure [2]) may indicate mass 
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motions with a very large velocity gradients, like shocks. During the flare, these lines exhibit 
broadened wings and red asymmetry, in addition to the line center shift. In addition, the Mg 
II subordinate lines are also significantly enhanced, providing more diagnostic power of the 
flaring atmosphere. IRIS shows a great potential to diagnose and help understand heating 
and dynamics of the lower atmosphere during solar flares. 

In the present study, we only focus on the spectra in the UV and EUV bands from two 
spectrographs, IRIS and EIS. In fact, for this XI.0 flare, there are also joint observations from 
some other instruments, such as the Solar Optical Telescope and the Dunn Solar Telescope, 
as well as AIA and the Helioseismic and Magnetic Imager. All these instruments provide 
rich data for us to study the different aspects of solar flares, which will be performed in a 
future work. 
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Table 1: Spectral Lines Used in This Stndy 


flDe Pontieu et al.l 120141 lYoung et al.l 120071) 


Ion 

Wavelength (A) 

log Tmax (K) 

Instrnment 

Mg II h 

2803.5 

4.0 

IRIS 

Mg II k 

2796.4 

4.0 

IRIS 

C II 

1334.5 

4.3 

IRIS 

C II 

1335.7 

4.3 

IRIS 

He II 

256.32 

4.7 

EIS 

Si IV 

1402.8 

4.8 

IRIS 

Fe XII 

192.39 

6.1 

EIS 

Fe XVI 

262.98 

6.4 

EIS 

Fe XXI 

1354.1 

7.0 

IRIS 

Fe XXIII 

263.76 

7.1 

EIS 
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Fig. 1.— Top panel: GOES 1-8 A soft X-ray light curve for the XI.0 flare. The grey and 
red vertical dashed lines represent the multiple IRIS and EIS scans over the flaring region, 
respectively. Bottom panels: IRIS slit-jaw 1400 A image (SJI 1400), AIA 1600 A image, and 
EIS He II 256 A intensity map several minutes after the GOES soft X-ray peak time. The 
red box shows the EIS FOV. 
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Fig. 2.— IRIS slit-jaw 1400 A images (column 1) and the slit spectra of Fe xxi, Si IV, C ll, 
and Mg ll h&k lines (columns 2-5) in the rise phase, at the peak time, and in the decay phase 
of the flare, respectively. In the slit-jaw images, the white lines indicate the slit positions, 
while the red bars show four ribbon pixels whose line prohles are plotted in Figures [3H6l 
respectively. The two green bars show the spatial range along the slit where the Doppler 
velocity and line width are presented in Figure [71 In the spectrum panels, the pink dotted 
lines mark the reference wavelength position of each line. 
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Fig. 3.— Line profiles of Fe xxi, Si IV, C ll, and Mg ll k at ribbon pixel NO marked in Figure 
121 The black curves refer to observed line profiles, and colored curves to fitted line profiles 
or flare-net/pre-flare emissions. The vertical dotted lines indicate the reference wavelength 
for each line. The Doppler velocities {V or 14) derived from Gaussian fitting are given in 
the Fe xxi and Si IV panels, and the velocities from bisector (14, at the ~40% maximum 
intensity marked by the diamond) or moment (14™) method given in the C ll and Mg ll 
panels. Negative and positive velocity values represent blueshifts and redshifts, respectively. 
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Fig. 4.— Line profiles of Fe xxi, Si iv, C ll, and Mg ll k at ribbon pixel SI marked in Figure 
[2J The black curves refer to observed line profiles, and colored curves to fitted line profiles 
or flare-net/pre-flare emissions. Note that there shows a saturation in the cores of the Si 
IV and C II lines (black dashed curves), and we analyze the unsaturated prohles closest to 
SI instead. The vertical dotted lines indicate the reference wavelength for each line. The 
Doppler velocities {V or V)) derived from Gaussian htting are given in the Fe xxi and Si iv 
panels, and the velocities from bisector ( 14 , at the ~40% maximum intensity marked by the 
diamond) or moment {Vm) method given in the C ll and Mg ll panels. Negative and positive 
velocity values represent blueshifts and redshifts, respectively. 
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Fig. 5.— Line profiles of Fe xxi, Si iv, C ll, and Mg ll k at ribbon pixel S2 marked in Figure 
[2J The black curves refer to observed line profiles, and colored curves to fitted line profiles 
or flare-net/pre-flare emissions. Note that there shows a saturation in the cores of the Si 
IV and C II lines (black dashed curves), and we analyze the unsaturated prohles closest to 
S2 instead. The vertical dotted lines indicate the reference wavelength for each line. The 
Doppler velocities {V or F)) derived from Gaussian htting are given in the Fe xxi and Si iv 
panels, and the velocities from bisector ( 14 , at the ~40% maximum intensity marked by the 
diamond) or moment {Vm) method given in the C ll and Mg ll panels. Negative and positive 
velocity values represent blueshifts and redshifts, respectively. 
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Fig. 6.— Line profiles of Fe xxi, Si iv, C ll, and Mg ll k at ribbon pixel S3 marked in Figure 
[21 The black curves refer to observed line profiles, and colored curves to fitted line profiles 
or flare-net/pre-flare emissions. The vertical dotted lines indicate the reference wavelength 
for each line. The Doppler velocities {V or 14) derived from Gaussian htting are given in 
the Fe xxi and Si IV panels, and the velocities from bisector (14, at the ~40% maximum 
intensity marked by the diamond) or moment (14n) method given in the C ll and Mg ll 
panels. Negative and positive velocity values represent blueshifts and redshifts, respectively. 
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Fig. 7 .— Evolution of Doppler (left columns) and non-thermal (right columns) velocities 
in the Fe xxi and Si iv lines along the slit indicated by the two green bars in Figure HI 
Overplotted on the Doppler velocity map in the Fe xxi line are the intensity contours of the 
Mg II h line, and those on the Si iv Doppler velocity map are the intensity contours of the Si 
IV line itself that is slightly saturated. The vertical solid line represents the GOES soft X-ray 
peak time. In the Doppler velocity panel of Fe xxi, the black data points indicate those 
Doppler blueshift velocities beyond the colorbar range. The box with dotted line shows the 
loop-top region that is used to obtain the average line center as the reference wavelength of 
Fe XXI. 
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Fig. 8.— Evolution of the intensity in four EIS spectral lines. The yellow line marks the 
GOES soft X-ray peak time. The time on the top row shows the starting time of each scan. 
All the intensity panels have the same FOV, as shown by the red box in Figure [H 
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Fig. 9.— Evolution of the Doppler velocity in four EIS spectral lines, superimposed on the 
intensity contours of the He ll line. The velocity range is ± 200 km s“^. The yellow line 
marks the GOES soft X-ray peak time. The time on the top row shows the starting time 
of each scan. The blank part in the Fe xxill and Fe xvi maps means no valnes, since the 
signal-to-noise ratio at those pixels is too low to get reliable htting resnlts. The dark pixels 
in the Fe xil velocity map represent a failed htting to derive the Doppler velocity dne to the 
Fe XII line snbmerged in the wing of the Fe xxiv line that becomes strong and broad aronnd 
the hare peak. 
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Fig. 10.— Sample line profiles (black curves) of four EIS emission lines observed during the 
flare, with the Gaussian fitting shown in colored curves. Note that, for the Fe xxill line in 
the top panels, we use either a single (right) or a double (left and middle) Gaussian function 
to £t the observed line prohles. The arrow marks a very weak blue asymmetry in the Fe 
XXIII line. The vertical dotted lines indicate the reference wavelength of Fe xxill. 
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Fig. 11.— Evolution of the Doppler velocity along the slit in four EIS lines. The slit position 
is around the same as the IRIS slit position shown in Figure [71 The yellow line marks the 
GOES soft X-ray peak time. The four labeled locations refer to places where the Doppler 
velocities are shown in Figure [TTl Note that, the Fe xxill velocities represent the shifts of 
blueshifted components, which are derived from a double Gaussian htting. The pixels within 
the red rectangular box in the Fe xxill panel show double-peak line prohles (see the top left 
panel in Figure [10]), and the other pixels only show a line prohle with a weak blue wing 
enhancement (see the top middle panel in Figure [TUD . 
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Fig. 12.— Doppler velocity as a function of temperature at four locations (marked in Figure 
[IID, derived from EIS EUV and IRIS UV lines. The dashed line refers to the zero velocity, 
and the two dotted lines show the range of velocity uncertainties. Negative velocity values 
represent blueshifts. Note that all the velocities are derived from a single Gaussian htting, 
except the two data points marked by a vertical bar in the Fe xxill line at locations 2 and 
3, which are derived from a double Gaussian htting. 














